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ABSTRACT 



A classification of ordinary differential equations and finite-difference 
equations in one variable having polynomial solutions (the gen- 
eralized Bochner problem) is given. The method used is based 
on the spectral problem for a polynomial element of the uni- 
versal enveloping algebra of s/2(R-) (for differential equations) or 
sl2(R)q (for finite-difference equations) in the "projectivized" rep- 
resentation possessing an invariant subspace. Connection to the 
recently-discovered quasi-exactly-solvable problems is discussed. 
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S. Bochner (1929) asked about a classification of differential equations 



Tip 



(0) 



where T is a linear differential operator of A;-th order in one real variable x G 
R and e is the spectral parameter, having an infinite sequence of orthogonal 
polynomial solutions (see 

Definition. Let us give the name of the generalized Bochner problem to 
the problem of classification of the differential equations (0) having (n + 1) 
eigenfunctions in the form of a polynomial of the order not higher than n. 

In Ref.2 it has been formulated a general method for generating eigenvalue 
problems for linear differential operators, linear matrix differential operators 
and linear finite-difference operators in one and several variables possessing 
the polynomial solutions. The method is based on considering the eigen- 
value problem for the representation of a polynomial element of the universal 
enveloping algebra of the Lie algebra in a finite-dimensional, 'projectivized' 
representation of this Lie algebra @. Below it is shown that this method pro- 
vides both necessary and sufficient conditions for the existance of polynomial 
solutions of linear differential equations and a certain class of finite-difference 
equations in one variable. 

1. Ordinary differential equations. 

Consider the space of all polynomials of order n 



where n is a non-negative integer and x G R. 

Definition. Let us name a linear differential operator of the A;-th order, 
Tfc quasi-exactly-solvable, if it preserves the space Vn- Correspondingly, 
the operator E^, which preserves the infinite flag Vq d Vi d V2 d ■ ■ ■ C 
P„ C . . . of spaces of all polynomials, is named exactly-solvable. 
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LEMMA 1. (i) Suppose n > [k — 1). Any quasi-exactly-solvable 
operator of k-th order Tk, can be represented by a k-th degree polynomial of 
the operators 

= x^dx — nx, 

j° = - 1 , (2) 

J~ = dx , 

(the operators (2) obey the sl2(R) commutation relations f\ ) . If n < {k — 1), 
the part of the quasi-exactly-solvable operator T^. containing derivatives up to 
the order n can be represented by a n-th degree polynomial in the generators 
(2). 

(a) Inversely, any polynomial in (2) is quasi-exactly solvable. 

(Hi) Among quasi-exactly-solvable operators there exist exactly-solvable 
operators Ek C T^. 

Comment 1. If we define the universal enveloping algebra Ug of a Lie 
algebra g as the algebra of all polynomials in generators, then at A; < n + 1 
is simply an element of the universal enveloping algebra f/si2(R) of ^^e algebra 
s/2(R) taken in representation (2). If > n + 1, then is represented as an 
element of f/si2(R.) P^^^ -B ^^^li , where B is any linear differential operator of 
the order not higher than {k — n — 1). 

Now let us proceed to the proof of the Lemma 1. 

Proof. 

Part I. One starts from the statement (i). Suppose n > (A; — 1). 

Step 1. It is easy to show that the coefficient functions standing before 
derivatives in any quasi-exactly-solvable operator should be polynomials. 

^The representation (2) is one of the 'projectivized' representations (see [^). 
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Step 2. Without loss of generality, any operator can be represented as 
a sum of homogeneous operators: 

n-j:o!^{x)di, (3) 

i=0 

where Ui indicates the degree of polynomial a^^^l'{x). Suppose, the maximal 
degree of homogenuity is M — ni — i > Now let us rewrite as sum of 
operators of fixed degree of homogenuity m: 

M k 

Tk-Y. ' Tt'' ^Y.^ti^'^^'dl (4) 

m=—k i=0 

(if m + z < 0, the corresponding A^^^ — 0). Evidently, if m < 0, then 
rpjjn) . •p^ •p^ Consider the operators Tj^"^^ at m > 0. Their action on 
monomials is 

T^""^ : ^ x'^"" (5) 
In order to preserve the space (1), the condition 

e + m< n , £ = 0,1,2, ...n (6) 

should be fulfilled. Generically, T^^"*'' is characterized by {k + 1) coefficients. 
The condition (4) leads to the system of m linear homogeneous equations for 
{k + 1) unknown coefficients ^4^™^ . Hence 

m<k (7) 

and the maximal degree of homogenuity of a quasi-exactly-solvable operator 
Tk is equal to k. 

Step 3. We will prove (i) inductively. 

(1). Take Tq. Evidently, Tq = const or, in other words, a zero-order 
polynomial in generators (2). 
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(2). Take Ti. As it follows from (4), (7) 

m=— 1 

We only need to consider the operator since the others preserve (1) 
automatically. In this case, the above-mentioned system of linear equations 
contain only one equation: 

A^il + nA^^l = , (9) 
which implies immediately, that 

ri'^=4;i^+- (10) 

Finally, 

t(^) = A^IJ^ + Aflj^ + + Afl + ^Afl (11) 

(^5). Now let us assume that a quasi-exactly-solvable operator Tfe_i is 
represented through generators (2). 

Take T^. In what follows from Step 2, the maximal degree of homogenuity 
of Tfc is equal to k. First consider Tj,''\ The corresponding system of linear 
equations on the coefficients A^/I^l.i = 0, 1,2, . . . A; contains k equations and 
hence has a non-trivial solution. This solution can be expressed through one 

(k) 

coefficient, e.g. A^ ^. It is easy to check that obtained coefficients correspond 
precisely to (J"'")'^ (cf.(9),(10)) and hence 

Tt^-4i{j+r (12) 

It is obvious, that (12) preserves the space (1). The remainder of is simply 
the A;-th order differential operator characterizing by degrees of homogenuity 
—k, —k+1, . . . ,0,1, . . . ,k — l. Now let us consider in the remainder the terms 



4 



containing the derivatives of the /c-order. Apparently, they can be rewritten 
in the following forms: 

^2k-3Qk ^ ^j+^k-2j0j- ^ ^^^2k-4Qk-l ^ _ _ _ 

9' = {J-f (13) 

(dots in (13) imply terms with lower order derivatives). The representation 
(13) is unambiguous, since there are only {2k + 1) different monomials of 
degree k in the generators J" (if we take into account the commuation re- 
lations and the elements of the ideal generated by the Casimir operator). 
Substituting (12), (13) into in the form (3), one gets 

fc+,fco,fc>0 

The operator T^-i can be represented through (2) by assumption. This ends 
the proof of the first part of (i) . 

Part II. Suppose n < k — 1. It is clear that the part of T^, containing 
derivatives of the orders (n + 1), (n + 2), . . . , /c annihilates the space (1) 
and these derivatives can keep any functions as coefficient functions. For the 
remainder of T^, containing derivatives up to the n-th order, the Part I of 
the proof holds. This concludes the proof of part (i) of the Lemma. 

Part III. The part (ii) of the Lemma is evident. The part (iii) of the 
Lemma is also evident from the proof of part (i) of the Lemma, since a quasi- 
exactly-solvablc operator, containing no operators of positive homogenuity, 
becomes exactly-solvable. □ 
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Comment 2. If the space (1) is considered not in general position 0, there exist 
hnear differential operators others than quasi-exactly-solvable ones, which 
map a certain polynomial of a degree n to some other polynomial of the 
same degree. 

Since 3/2(1^) is a graded algebra, let us introduce the grading of generators 

(2): 

deg{J+) = +1 , rfec/(J°) = , deg{J-) = -1, (15) 

hence 

(ie^[(J+)"+(J°)"«(J-)"-] = n+-n^. (16) 
The grading allows to classify the operators Tk in Lie-algebraic sense. 

LEMMA 2. A quasi-exactly-solvable operator Tk C ^/^ZjCR) ^^■^ 
terms of positive grading, iff it is an exactly-solvable operator. 

Comment 3. It is easy to see that the grading is nothing but the homogenuity, 
which has been introduced in the proof of the Lemma 1 and the statement 
of this Lemma becomes obvious. 

Definition. Let us call the operator T{x) symmetric, if one can introduce 
the scalar product 

{f,9)p = / f{x)g{x)p{x)dx 
Jr 

for positive p > 0,p{x) G S{R), such that 

iTf,g)p = {f,Tg), 
where the Schwatz space is defined as 

S{R) = {p:R^R,\ p('-)(x) |< Ckn{l+ I X I)-", VA;, n} . 

^ 'not in general position' means that, it is allowed to have some correlations between 
coefficients of original polynomial and an operator Tk- Otherwords, not any polynomial 
of degree not higher than n is mapped to Vn- 
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THEOREM 1. Let n is non-negative integer. Take the eigenvalue problem 
for a linear differential operator of the k-th order in one variable 

Tkif = eip , (17) 

where Tk is symmetric. The problem (17) has {n + 1) linear independent 
eigenfunctions in the form of polynomial in variable x of the order not higher 
than n, iff Tk is quasi-exactly-solvable. The problem (17) has an infinite 
sequence of polynomial eigenfunctions, iff the operator is exactly-solvable. 

Comment 4- The "if part of the first statement is obvious. The "only if 
part is a direct corollary of Lemma 1. 

This theorem gives a general classification of differential equations 

k 

Y,ajix)ip^^\x) = eip{x) (18) 

j=0 

having at least one polynomial solution in x, thus resolving the generalized 
Bochner problem. The coefficient functions aj{x) must have the form 

k+j 

= (19) 

The explicit expressions (19) for coefficient function in (18) are obtained by 
the substitution (2) into a general, A;-th degree polynomial element of the 
universal enveloping algebra Ugi^^ny Thus the coefficients aj^i can be ex- 
pressed through the coefficients of the A;-th degree polynomial element of the 
universal enveloping algebra t/s/2{R)- The number of free parameters of the 
polynomial solutions is defined by the number of parameters characterizing 
a general k-th degree polynomial element of the universal enveloping algebra 
Usi2(R) 0- Rather straightforward calculation leads to the following formula 

par{Tk) = {k + If (20) 

Counting free parameters, one should introduce a certain ordering of generators to 
avoid double counting because of commutation relations. Also the quadratic Casimir 
operator and the double-sided ideal generated by it should not be taken into account. 
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where we denoted the number of free parameters of operator Tk by the symbol 
par{Tk). For the case of an infinite sequence of polynomial solutions the 
expression (19) simplifies to 

j 

= Y.^i,i^' (21) 

in agreement to the results of H.L.Krall's classification theorem 0(see also 
In this case the number of free parameters is equal to 

In present approach Krall's theorem is simply a description of differential 
operators of A;— th order in one variable preserving a finite flag Vq G Vi G 
V2 G . . . G Vk oi spaces of polynomials. One can easily show that the 
preservation of such a set of polynomial spaces implies the preservation of an 
infinite flag of such spaces. 

One may ask a more general question: which non-degenerate linear dif- 
ferential operators have finite-dimensional invariant sub-space of the form 

. . . , , (23) 

where a{x) is a function and I3{x) is a diffeomorphism of the line. Such 
operators are obtained from the ones of the Theorem 1 by the change of 
variable x 1— > (3{x) and the "gauge" transformation ip{x) ^ a{x)tfi{x)). 

Let us consider the set of second order differential equations (17), which 
can possesses polynomial solutions. From Theorem 1 it follows that the 
corresponding differential operator must be quasi-exactly-solvable and can 
be represented as 

T2 = C++ J+ J+ + c+o J+ J° + C+-J+J- + co-J^r + C—J-J- + 

c+J^ + CqJ^ + c^J- + c, (24) 
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where Ca/s, Ca, c e R. The number of free parameters is par{T2) — 9. Under 
the condition C++ = c+o — c+ — 0, the operator T2 becomes exactly-solvable 
(see Lemma 2) and the number of free parameters is par{E2) — 6. 

LEMMA 3. // the operator (21) is such that 

Tl 

C++ = and c+ = (— — m)c+o j m = 0, 1, 2, . . . (25) 

then the operator T2 preserves both Vn and Vm- The number of free param- 
eters is par{T2) = 7. 

PROOF. By straightforward analysis. 

In fact, the Lemma 3 means that T'2( J°(n), Cq,^, c^) can be rewritten as 
J"(m), c^^, c^). As a consequence of Lemma 3 and Theorem 1, among 
polynomial solutions of (17) there are polynomials of order n and order m. 

Remcirk. Prom the Lie-algebraic point of view Lemma 3 means the ex- 
istance of representations of second-degree polynomials in the generators (2) 
possessing two invariant sub-spaces. In general, if n is not a non-negative in- 
teger in (2) ( correspondingly, (2) becomes infinite-dimensional), then among 
representations of A;-th degree polynomials in the generators (2), lying in the 
universal enveloping algebra, there exist representations possessing 0, 1, 2, /c 
invariant sub-spaces. Also this properly implies existence of representations 
of the polynomial elements of the universal enveloping algebra, which can be 
obtained starting from different representations of the original algebra. Even 
starting from an infinite-dimensional representation of the original algebra, 
one can construct the elements of the universal enveloping algebra having 
finite-dimensional representation (e.g. the parameter n in (25) is non-integer, 
however T2 has the invariant sub-space of the dimension {m-\- 1)). 

Substituting (2) into (24) and then into (17), we obtain 

P^{x)dl^{x) + P^{x)d,^{x) + P2{x)^{x) = e^{x), (26) 
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where the Pj{x) are polynomials of j'-th order with coefficients related to 
Ca/3,Ca and n (see (9)). In general, problem (26) has {n + 1) polynomial 
solutions. If n = 1, as a consequence of Lemma 1, a more general spectral 
problem than (26) arises 

F3{x)dlip{x) + Q2{x)dMx) + Qi{x)ip{x) = eip{x), (27) 

possessing only two polynomial solutions of the form {ax + b), where is 
an arbitrary real function of x and Qj{x),j = 1,2 are polynomials of order 
j. For the case n = (one polynomial solution) the spectral problem (11) 
becomes 

F2{x)dl^{x) + Fi{x)dMx) + Qofix) = e^{x), (28) 

where F2^i{x) are arbitrary real functions of real x and Qo is a real constant. 
After the transformation 

t : ip ^ ^{x{z))e^^'\ (29) 

where z h-* x{z) is a diffeomorphism of the line and A{z) is a certain real 
function, one can reduce (26)-(28) to the Sturm-Liouville problem 

{d! + Viz))ip = Sif, (30) 

with the potential 

V{z) = {A'f + A" + P2{x{z)) 

where A = J(^)dx~logz' for (23). If the functions (29), obtained after trans- 
formation, belong to the L2(I')-spaceQ , we reproduce the recently discovered 
quasi-exactly-solvable problems @] , where a finite number of eigenstates was 
found algebraically. For example, 

T2 = -4J°J- + 4a J+ + 46 J° - 2{n + 1)7" (31) 

^In dependance on a diffeomorphism z x{z), the space V can be the infinite real 
Une, semi-infinite real line and a finite real interval. 
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leads to the spectral problem (30) with the potential 

V{z) = a^z^ + 2abz^ + {b^ - {An + 3)a)z^, (32) 

for which the first (n+1) eigenfunctions, even in x can be found algebraically. 

It is worth noting that the use of (27) as input leads to one-functional 
family of the Schroedinger operators with two explicitly known eigenstates. 
One of such operators has been described at 0, where the authors confus- 
ingly stated about non-existance of s/2(R.) algebra behind. 

Taking different exactly-solvable operators E2 for the eigenvalue problem 
(17) one can reproduce the equations having the Hermite, Laguerre, Legendre 
and Jacobi polynomials as solutions 0. Also under special choices of 
general element i?^, one can reproduce all known fourth order differential 
equations giving rise infinite sequences of orthogonal polynomials (see 
and other papers in this volume). 

Recently, A.Gonzalez-Lopez, N.Kamran and P.Olver @] gave the com- 
plete description of second-order polynomial elements of t/s/2{R) leading to 
the square-integrable eigenfunctions of the Sturm-Liouville problem (30) af- 
ter transformation (29). Consequently, for second-order ordinary differential 
equations (26) the combination of this result and Theorem 1 gives a gen- 
eral solution of the Bochner problem as well as the more general problem of 
classification of equations possessing finite number of orthogonal polynomial 
solutions. 

2. Finite-difference equations in one variable. 

The generalized Bochner problem is defined in the same way, as it has 

been done for differential equations. The only difference is to consider the 

operator T in the problem (0) as a linear finite-difference operator of a finite 

^For instance, putting the parameter a = in (31), the equation (26) converts to the 
Hermite equation (after some substitution.) 
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order. For the case of one real variable a solution of the classification problem 
is very similar to the case of ordinary differential equations described above. 

Let us introduce the finite-difference analogue of the operators (2) ^ 

= x^D — {n]x 

J0= xD-h (33) 
J' = D, 

where h = '^"^^""^^-^ , {n} = is the quantum symbol, q is a number 

characterizing the deformation, Dz = 1 + qzD and Df{z) = '^''^(IZ^-jl^^ + 
f{qz)D is a shift or a finite-difference operator (or so called the Jackson 
symbol (see 0)). The operators (33) after multiplication by some factors as 

~ ^ g-" {2n + 2} J, 
p+1 {n+1} 

( see 0) form a quantum s/2(R)g algebra with the following commutation 
relations 

ifr - j~f = -j' 
q'j^r - rj+ = -(g + i)? (34) 

?j+ - qrr = r 

(this algebra corresponds to the second Witten's quantum deformation of s/2 
in the classification of CZachos [^]). If g ^ 1, the commutation relations (34) 
reduce to the standard s/2(R.) ones. A remarkable property of generators (33) 
is such that, if n is a non- negative integer, they form the finite-dimensional 
representation. 

Similarly as for differential operators one can introduce quasi-exactly- 
solvable and exactly-solvable operators. 
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Definition. Let us name a linear difference operator of the k-th order, 
quasi-exactly-solvable, if it preserves the space Vn- Correspondingly, the 
operator Ek, which preserves the infinite flag Vo C Vi C V2 C . . . C Vn ■ ■ ■ 
of spaces of all polynomials, is named exactly-solvable. 

The analogue of the Lemma 1 holds. 

LEMMA 4. (i) Suppose n > {k — 1). Any quasi-exactly-solvable 
operator of k-th order Tk, can he represented by a k-th degree polynomial 
of the operators (33). If n < {k — 1), the part of the quasi-exactly-solvable 
operator containing derivatives up to the order n can be represented by a 
n-th degree polynomial in the generators (33). 

(a) Inversely, any polynomial in (33) is quasi-exactly solvable. 

(Hi) Among quasi-exactly-solvable operators there exist exactly- solvable 
operators Ek dT^. 

PROOF. Straightforward analogue of the proof of the Lemma 1. 

Comment 5. If we define an analogue of the universal enveloping algebra 
Ug of a Lie algebra g as an algebra of all polynomials in generators sl2(R)q. 
Then a quasi-exactly-solvable operator at A; < n -|- 1 is simply an element 
of the 'universal enveloping algebra' Usi^cr)^ of the algebra sl2(R)g taken in 
representation (33). If A; > n -\- 1, then is represented as an element of 
Usi2(R)g plus BD^'^^, where B is any linear difference operator of the order 
not higher than {k — n — 1). 

Similar to s/2(R) , one can introduce the grading of generators (31) 
of s/2(R)g (see (15)) and, hence,of monomials of the universal enveloping 
Usi2{R)q (see (16)). The analogue of Lemma 2 holds. 

LEMMA 5. A quasi- exactly- solvable operator C Usi2{-R)q has no 
terms of positive grading, iff it is an exactly-solvable operator. 
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PROOF. Straightforward analogue of the proof of the Lemma 2. 

THEOREM 2. Let n is non-negative integer. Take the eigenvalue prob- 
lem for a linear difference operator of the k-th order in one variable 

f^^{x) = eip{x), (35) 

where Tk is symmetric. The problem (35) has {n + 1) linear independent 
eigenfunctions in the form of polynomial in variable x of the order not higher 
than n, iff is quasi-exactly-solvable. The problem (35) has an infinite 
sequence of polynomial eigenfunctions, iff the operator is exactly-solvable E^. 

Comment 6. Saying the operator Tj^ is symmetric, we imply, that consid- 
ering action of this operator on a space of polynomials of degree not higher 
than n, one can introduce a positively-defined scalar product and the oper- 
ator Tfc is symmetric with respect to it . 

Comment 7. The "if part of the first statement is obvious. The "only 
if part is a direct corollary of Lemma 4. 

This theorem gives a general classification of finite-difference equations 

k 

Y,aj{x)D^ ip{x) = eip{x) (36) 

having polynomial solutions in x. The coefficient functions must have the 
form 

k+j 

o-A^) = Y.V'- (37) 

Particularly, this form occurs after substitution (33) into a general k-th. de- 
gree polynomial element of the universal enveloping algebra Usi2{ii)q- It guar- 
antees the existance of at least a finite number of polynomial solutions. The 
coefficients related to the coefficients of the k-ih. degree polynomial 
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element of the universal enveloping algebra Usi2(R)g- The number of free pa- 
rameters of the polynomial solutions is defined by the number of free param- 
eters of a general k-th order polynomial element of the universal enveloping 
algebra Usi^i^r)^^ A rather straightforward calculation leads to the following 
formula 

par{fk) = {k + lf + l 

(for the second order finite-difference equation par{T'^) = 10). For the case 
of an infinite sequence of polynomial solutions the formula (37) simplifies 

j 

dj{x) = ^^cLj^iX^ (38) 

1=0 

and the number of free parameters is given by 

(k + l)(k + 2) 
par{Ek) = ^ 'j + 1 

(for k = 2, par{E'^) = 7). The increase in the number of free parameters 
compared to ordinary differential equations is due to the presence of the de- 
formation parameter q. In [§] one can find the description in present approach 
of the g-deformed Hermite, Laguerre, Legendre and Jacobi polynomials (for 
definition these polynomials see 0). 

The analogue of Lemma 3 holds as well 

LEMMA 6. // the operator Ta (see (24)) is such that 

C++ = and c+ = {n — {m})c+o , at m = 0,1,2, .. . (39) 

^For quantum sZ2(R-)i3 algebra there are no polynomial Casimir operators However, 
in the representation (33) the relationship between generators analogous to the quadratic 
Casimir operator 

qJ+J- - J° J° + {{n + 1} - h)j" = h{n - {n + 1}) 

appears. It reduces the number of independent parameters of the second-order polynomial 
element of C/si2(R.)q- becomes the standard Casimir operator at p ^ 1. 
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then the operator T2 preserves both P„ and Vm and polynomial solutions in 
X with 8 free parameters occur. 

PROOF. By straightforward analysis. 

Rather interesting situation occurs, if the parameter of deformation q 
equals to the root of unity. 

LEMMA 7. // a quasi-exactly-solvable operator preserves the poly- 
nomial space Vn and the parameter q is satisfied to the equation 

= 1 , (40) 

then the operator preserves an infinite flag of polynomial spaces Vq C 

C T'sn C. . . C Pfc„ C . . .. 

It is worth emphasizing, that in the limit q stresses to one. Lemmas 4,5,6 
and Theorem 2 coincide to the Lemmas 1,2,3 and Theorem 1, respectively. 
Thus the case of differential equations in one variable can be considered 
as particular case of finite-difference ones. Evidently, one can consider the 
finite-difference operators, which are a mixture of generators (33) with the 
same value of n and different g's. 

In closing, I would like to thank to L.Michel, V.Ovsienko, S.Tabachnikov 
and, especially, V.Arnold and M.Shubin for numerous useful discussions. 
Also I am very grateful to the Centre de Physique Theorique for hospitality 
extended to me, where this work has been completed. 
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